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ANALYSIS OF TEMPERATURE DISTRIBUTION IN LIQUID-COOLED TURBINE BLADES 1 

By John X. B. Livincood and W. Btron Brown 


SUMMARY 

The temperature distribution in liquid-cooled turbine blades 
determines the amount of cooling required to reduce the blade 
temperature to permissible values at specified locations. This 
report presents analytical methods for computing temperature 
distributions in liquid-cooled turbine blades, or in simplified 
shapes used to approximate sections of the blades. 

The individual analyses are first presented in terms of their 
mathematical development. By means of numerical examples, 
comparisons are made between simplified and more complete 
solutions and the effects of several variables are examined. 
Xondimensional charts . to simplify some temperature- 
distribution calculations are also given. 

It was found that for blade materials having low thermal 
conductivities ( corresponding to blade materials in current use) 
a one-dimensional s panwise temperature-distribution equation 
is applicable near a coolant passage. The cooled part of the 
blade, irrespective of blade length, was found to remain at a 
determinable uniform temperature , called the prevalent blade 
temperature. The prevalent blade temperature increased about 
160° F for an increase in the difference between the effective gas 
temperature and the coolant temperature of 1000° F for a range 
of effective gas temperature from 2000 ° to 6000° F. 

In cases where rim cooling is insignificant, it was found that a 
one-dimensional chordwise temperature-distribution equation 
was sufficiently near the relaxation solution for the actual blade 
shape for a first-order approximation. 

INTRODUCTION 

A limitation on design and performance of aircraft gas- 
turbine power plants is the strength of the turbine materials, 
which decreases as temperature increases. This limitation 
may be greatly alleviated, even when nonstrategic materials 
are used, by the application of turbine cooling. An extensive 
study of both liquid and air cooling of turbine blades has been 
carried out at the NACA Lewis laboratory since 1945. 

The temperature distribution in a turbine blade determines 
the amount of cooling required to reduce the blade tempera- 
ture to permissible values at specified locations. Analytical 
studies have been conducted to predict temperature distribu- 
tions in cooled turbine blades. A summary of the analytical 
methods developed for calculating spanwise temperature 
distributions for three types of air-cooled turbine blade is 
presented in reference 1. The present report summarizes 
the analytical methods developed for calculating temperature 
distributions in liquid-cooled turbine blades. Parts of these 
analyses are presented in references -2 and 3. 


Because the calculation of a generalized three-dimensional 
temperature distribution through a liquid-cooled turbine 
blade is too tedious to be of practical value and because some 
knowledge of the distribution in a radial direction is neces- 
sary, the first investigation presented herein is limited to a 
one-dimensional spanwise distribution. This distribution, 
however, is valid only in the immediate vicinity of the 
coolant passages. In order to study more accurately the 
temperature distribution in the trailing section of the blade, 
a three-dimensional solution is determined for a rectangular 
parallelepiped used to approximate a blade trailing section. 
This solution serves to determine the importance of blade 
thermal conductivity and of distance from the coolant pas- 
sage. Because numerical calculations for a rectangular 
parallelepiped produce a constant spanwise temperature in 
the region of the blade beyond the influence of rim cooling 
(that is, in the region of the blade unaffected by conduction 
to the rim), one-dimensional chordwise temperature distribu- 
tions are then calculated for shapes approximating various 
parts of a blade cross section; these investigations are valid 
for the part of the blade beyond the influence of rim cooling. 
Finally, a two-dimensional temperature distribution through 
the cross section of an actual blade shape is obtained in order 
to determine the accuracy of the previously determined 
approximate solutions and to study the effect on temperature 
distribution of a peripheral variation of gas-to-blade heat- 
transfer coefficient. 

Numerical examples based on coefficients available from 
unclassified sources are included. 

METHODS OF ANALYSIS 

One-dimensional spanwise, three-dimensional, one-dimen- 
sional chordwise, and two-dimensional chordwise temperature- 
distribution equations for liquid-cooled turbine blades 
or for simplified shapes used to approximate sections of 
liquid-cooled turbine blades are presented. For simplicity 
in presentation, the individual analyses are first formulated 
in terms of the mathematical methods. The comparative 
applicabilities of the methods of analysis are established in 
later sections of the report. 

For all the analyses presented, the following conditions are 
assumed : 

(1) The coolant temperature is constant at the average 
temperature of the coolant, or the coolant forms a constant- 
temperature reservoir. This assumption is valid when the 
change in coolant temperature is small in comparison with 
the temperature difference between the gas and the coolant. 


1 Supersedes NACA TN 2323, “Analysis ot Temperature Distribution In Liquid-Cooled Turbine Blades” by John N. B. Livfngood and W. Byron Brown, 1951. 
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(2) The blade-to-coolant heat-transfer coefficient is con- 
stant. 

(3) The thermal conductivity of the blade is constant. 

(4) The effective gas temperature at all parts of the blade 
is the same. 

(5) Radiation effects are considered to be included in. the 
lleat-transfer coefficients. The following individual analyses 
are presented : 

One-dimensional spanwise temperature distribution, — A 
one-dimensional spanwise calculation gives a radial tempera- 
ture distribution valid in a part of the blade near a coolant 
passage. In this particular case, the radial temperature 
distribution is carried through the blade and tlirough the 
turbine rotor. Heat transfer, between the turbine rotor and 
the fluid on the outside of the rotor is assumed to take place 
and a constant value of gas-to-blade heat-transfer coefficient 
is assumed. For simplicity, a blade of constant cross- 
sectional area and perimeter is considered. 

For blades with little taper, the method is applicable if 
average cross-sectional area and perimeter are considered. 
In general, however, a numerical solution is necessary for 
tapered blades. 

Three-dimensional temperature distribution through sim- 
plified trailing section. — A. three-dimensional temperature 
distribution for the trailing section of the blade is obtained 
by considering this part of the blade as a rectangular 
parallelepiped. The blade-root temperature is considered 
constant, the gas-to-blade heat-transfer coefficient is constant, 
and it is assumed that no heat flows across the median plane 
of the simplified blade section. 

An approximate solution for tapered blades is obtainable 
by this method if the blade is considered in small sections, 
average values of cross-sectional area and perimeter are 
considered for each section, and appropriate compatibility 
relations are satisfied. 

One-dimensional chordwise temperature distribution 
through simplified shapes. — Rectangular and trapezoidal 
sections, concentric circle annuli, and sections between 
parallel plates are considered as the simplified shapes used 
to approximate parts of a turbine-blade cross section. 
The gas-to-blade heat-transfer coefficient is considered 
constant. For comparison purposes, the rectangular and 
trapezoidal sections are constructed so that lengths and 
areas are identical. 

Two-dimensional chordwise temperature distribution, — 
Numerical two-dimensional temperature distributions 
through the cross section of an actual blade shape are de- 
termined. Constant and variable gas-to-blade lieat-transfer 
coefficients are considered. 

ONE-DIMENSIONAL SPANWISE TEMPEBATUBE DISTBIBUTION 

The spanwise temperature-distribution equations for a 
liquid-cooled turbine, a section of which is shown in figure 1, 
are derived in reference 2 and are reproduced in appen- 
dix B. In this spanwise case, the investigation was extended 
from blade tip to rotor hub. 

For convenience, the turbine was divided into four 
sections (fig. 1), and the temperature-distribution equations 
were obtained from heat balances for differential elements 



Uncooled section of 
blade (section I) 

Cooled section of 

blade (sections) 

Rim (section 3) 


Rotor (section 4) 


. Figure 1.— Arrangement of Internal-cooling passages In turbine. 


in each section. Inasmuch as chordwise conduction w'as 
neglected, the following equations, valid only in the neighbor- 
hood of the coolant passages, were obtained. For the 
uncooled .section of the blade (fig. 1, section 1), it was found 
that 

— Ci cosh A(xi+C t ) (1) 


■where C x and C 3 are integration constants and 



(All symbols are defined in appendix A.) The value of C 4 
is determined by use of the boundary condition at the blade 
tip (where heat enters by convection only). 


or 


— IcbAb, i — h 0 A B , i {T ,, t — T b , t ) 


tanh AC 2 = 


h/Q 

O 


V U : s / 


For the cooled section of the blade (fig. 1, section 2), it 
was found that 


T B ,i=T Bi p—Cze Br 2—Cie Bx i 

where C t and C\ are integration constants, 

rp iTl 

hJo+h t J t . a 

and 

\ k B A B . 2 ) 


( 2 ) 


(2a) 
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The prevalent blade temperature T B , r is the temperature 
the blade would assume if no heat flows in or out of the blade 
ends. 

For the rim section (fig. 1, section 3), with assumed con- 
stant area, 

T^G+C^+C^h (3) 

where C s and are integration constants, 

4 fl~r s,aihgT a -{-(hi'J.t'i)a,Tr \ 

gj |~ ^^3, a r ^q ~t~ t l)qg ~l i 

L 1" rlm^lrlm J 

and T a is the fluid temperature on the outside of the rotor 
and (h,i,ilt,t)<u is the average of the values of found 

for the blade coolant passages, the rim circumferential 
passages, and the radial passages through the rotor (fig. 1). 

For the rotor section ffig. 1, section 4), with assumed con- 
stant strength, 

Tr=J -f C- cosh H r (4i 

where C 7 is an integration constant, and where 

j_ 4irr i.aihxTg-^iht^li^ apTi 
■Ivl* 4,as^ad~( at 

and 

rr ["~ 4 IT T Ai.. i)ai~ ]1 

L k B A B J 

\'alues of the six intergration constants C u C' t , C ir C\, C e , 
and C 7 are found by solving simultaneously the six equations 
resulting from equating temperatures and heat flows at the 
various junction points of the four sections of the turbine. 
(Subscripts a and |3 used after numerical subscripts designate 
the end nearest the blade tip and the end nearest the rotor 
hub of the sections to which the numerical subscripts refer, 
respectively.) At the junctions of sections 1 and 2, 

T t , t - C\ cosh A{x u p+ CJ = T B ,,-C 3 e Bx > °- C,e~ B A« ( 5 ) 

k B A BilJ CiA sinh A(x u( i+ C^)=Bk B A Bw2r a (C t e Bz *- “ — 

C* ~ BXi - -) +htJA B ' A b , , a ) (T b , , - C,e B A « - 

r 4 e-^--r,) (6) 

Equation (6) equates the heat leaving section 1 to the sum 
of the heat entering the metal of section 2 and the heat 
entering the coolant at that part of the blade where the 
inlet and outlet passages are connected. An approximation 
in the procedure has been introduced at this point because 
a separate blade section for the part of the blade containing 
the connecting passage between inlet and outlet passages is 
not introduced. A numerical calculation showed that use 
of such a section would slightly decrease the temperature at 
the blade tip. 


At the junction of sections 2 and 3, 

' T b , , - C 3 e B ^- ’ - C,e~ Bz ^= G + C,e E A « + <? 6 « _£r * « (7) 

k B A Bti B(C 3 e Bl *->- CVT B ^)+ 

A^-ZA^r h 0 A Bti (T t ,.- T B .,+ C t e B '*'+ 

Zj j&b, 5 

C^e - Bl *- s ) = Ek B 2 (C’ 8 «" £f3 '“-C , 5 e &s “) (8) 

Equation (8) equates the sum of the heat leaving section 
2 and that entering section 3 directly from the hot gases to 
the total heat entering section 3. 

At the junction of sections 3 and 4, 

G J rC i e Els '<‘+C a e~ Elt -i , =J +C 7 cosh Hr a (9) 

k^Ei-C^o +C a e~^)=4^- k B HC 7 sinh Hr a (1 0) 

^-rlm 

THREE-DIMENSIONAL TEMPERATURE DISTRIBUTION THROUGH 
SIMPLIFIED TRAILING SECTION 

The spanwise temperature-distribution equations previ- 
ously presented are valid only in the neighborhood of the 
coolant passages. Most currently used turbine blades are 
so shaped that coolant passages cannot be located near the 
trailing edge. Because of the impossibility of placing reason- 
ably large coolant passages near the trailing edges of con- 
ventional turbine blades and because the trailing sections 
seem most likely to be the hottest parts of the blade, other 
detailed studies were made of temperature distributions. 
The first study was devoted to approximating the trailing 
section by a rectangular parallelepiped and to determining 
the three-dimensional temperature distribution through this 
parallelepiped (fig. 2). 

(= j i 



Fjgusx 2.— Rectangular parallelepiped used to approximate trailing section of 4-Ineh turbine 
blade for three-dimensional analysis. 
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The median plane of the rectangular parallelepiped was 
chosen as the plane 3=0, and symmetry was assumed about 
this plane (no heat flow was assumed to cross this plane, 
as previously stated). Boundary conditions at the blade 
trailing edge and blade tip were simplified by assuming the 
blade width and blade length to be extended by a distance 
equal to one half the blade thickness. In figure 3, the edge 
AIN at temperature T gained some heat from the gas stream. 
The extended surfaces NO and MP were at nearly the same 
temperature T and no heat entered the edge OP; these 
surfaces therefore gained practically the same amount of 
heat as the actual exposed edge. The validity of this ap- 


Liquid 
coolant 
at T, 


'//A 


< 

r 

j 



j 

1 i 

M P 


.... i 




N 0 

Figure 3.— Correction for heat reoeived by trailing edge. 


Values of K r 


and O m „ are given by the relations 
K 


(-1)- 1 sm Nm 2 8h< (: T t . t - Ti) 


\(2n — 1) 


and 


proximation is discussed in reference 4 (pp. 216 and 217). 
Distances increased by r/2 are denoted by primes. 

The derivation of the three-dimensional temperature dis- 
tribution, originally derived in reference 3, is reproduced in 
appendix C. The differential equation, in final form, for 
the temperature distribution is found to be 


b*d , b 2 9 , b 2 9 
bx' 2+ by' 2+ b3 2 


= 0 


(ID 


where x ' and y' denote x and y increased by r/2, respectively. 
A solution of equation (11) satisfying the boundary condi- 
tions at x'—y'=z=0 is 


B=±, £(«».» cos L n x' cosh M n _ „ y' cos N m z + 

m — 1 t »— 1 

O m , n cosh P M ,»x' cos Q„y' cos N m z) 


( 12 ) 


where if, L, M, N, 0, i P, and Q are constants. Relations 
among the constants are 


M mta *=NJ+LJ 

P*. n 2 =N m *+Q 


t \m,n= 

n ) 


1,2,3, 


(13) 


N m \ I sin 2N m 


( 


r\ \k B M m-n sinh M m , n j'+hi cosh M 


m,nj ) 


(14) 


sin 


O m . 


sin Q n j’ AjT^.-T^r) 

Qnj' (, | sin 2Q n j' \ N t 
V + 2 0 n j' ) " n 2 


/ sin 2iy,| \ 

l / 


cosh P n , n V 


(15) 


Application of the boundary conditions at x'=b', y'=j', and 2 =r/ 2 lead to the following methods of evaluating L, Q, and 
N, respectively: 

L *=(n- 4)f (16) 


, /% •/ h{ kg 

tan H ,3 


( 17 ). 


and 


h„T 


tan N m 


t _ h 0 2 k B 

2 k B N n N n r 


-2 



It is at once obvious that from equation (16) many values 
of L n result. Because of the periodicity of the tangent 
function, equations (17) and (18) likewise have many 
solutions. 

Values for all the constants in equation (12) can thus be 
found and the temperature can now be computed at any 
point in the rectangular parallelepiped. In reference 3, 
appendix B, it is shown that sufficiently accurate results can 
be obtained by using m = 1 only. 


ONE-DIMENSIONAL CHOBDWISE TEMPERATURE DISTRIBUTIONS 
THROUGH SIMPLIFIED SHAPES 

Because the spanwise temperature distributions are valid 
only in the neighborhood of the coolant passages and because 
the three-dimensional approximate solution resulted in a 
constant spanwise blade temperature in the part of the blade 
beyond the influence of rim cooling, one-dimensional chord- 
wise temperature distributions were determined for sections 
of a liquid-cooled turbine blade that can be approximated 
by simple shapes (fig. 4). Rectangular and trapezoidal 
approximations for blade trailing sections were considered. 
Trapezoidal sections may also be used to approximate the 
leading section in some blades. In addition, analyses were 
made for leading sections approximated by concentric circle 
annuli and for the sections of blades with very little metal be- 
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(a) Rectangle. 

(b) Trapezoid. 

(Cl Concentric circles. 

(d) ParaUel plates. 

Ftot Hl 4.— Simplified shapes used In ono-dimenstonal chordwlse analysis. 


t ween blade outer surface and coolant passages approximated 
as regions between parallel plates. The temperature distri- 
butions through, these simplified shapes, derived in appendix D 
and valid in regions beyond the influence of rim cooling, follow. 

Rectangular trailing section of blade removed from 
influence of rim cooling. — For a blade trailing section 
approximated by a rectangle (fig. 4(a)) the temperature- 
distribution equation is 


where 

and 


6=T,.'-T B 
( 2 k 0 \i 


Trapezoidal trailing section of blade removed from in- 
fluence of rim cooling. — For a trapezoidal approximation to 
the blade trailing section (fig. 4(b)), the temperature distri- 
bution is 


jSk B iTt ’ ~ Tl) J ° {i « +* f)l 
[ i Ji (» f o ff ,(i r $■ - [h ,(i r o i Ji (» r *>i ■ + dSt r 


( 20 ) 


where 


( h y 

Vi's sin i pj 


K 


ip=tan _1 


Tj~T l 

2 J 


r=Zfi(ifi)«7..(if»)+t«7i(iri)iH tf (ira) 

and J e , ij u H 0 , and iH x are Bessel functions. For a wedge- 
shaped section the temperature-distribution equation (20) 
reduces to 

Mb 

. (21) 

Concentric-circle-annulus approximation for blade leading 
section. — The equations expressing the blade-metal tem- 
perature at the coolant passage wall T Bf t, at the leading edge 
of the blade T B , 0I and at any point in the annulus T B are 
found to be (fig. 4(c)) 

2 k B 

kzv.-r,.,) 


I 


and 


2k B 

Mt)J 


T b .o — 7f.« — < 


D 0 h 


( 22 ) 

W\ IK T B ,,-T B , t ) (23) 


i0ga (^) 

r rr I \*Si/ (rp rp \ 

B — 1 B,i ‘ 1 ,-r\ \ UB,o — J-BA) 

l0 *'(w 


-p (T,, t —Ti) cosh 
d- B 


where 

t b Tr ,= r ^ t * r ^‘ 

■19) 

B ‘° ±B - t , , t B f 2 , 2 \ 

q sinh <pj'+jr cosh <pj' 



(24) 

(25) 


Section of blade approximated by parallel plates. — The 
equations expressing the blade-metal temperatures at the 
coolant-passage wall (fig. 4 (d)) T B%i and at the blade outei 
edge T b , 0 are 

r,+(^>.n.+w 

1+ Gxx) (i * +A,> 


(26) 


and 


T b .o=- 


where o- is the distance between the parallel plates. 


(21 
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TWO-DIMENSIONAL CHOBDWISE TEMPERATURE DISTRIBUTION 


Tlie temperature at any point on the blade cross section 
may be found by solving the Laplace differential equation 
given in terms of the temperature difference 6 


b 2 8 , b*d 
dX 2+ dT- 


0 


(28) 


where A” and Y are the Cartesian coordinates in the blade 
cross section. The boundary conditions, expressed in 
terms of the partial derivative of 6 in the direction normal 
to the boundary, are 


do h„ 
dn F b 


(29) 


at the outer boundary and 


d0 _hj 
dn~k B 




(30) 


at the coolant-passage boundary'. 

A closed solution to equation (28) cannot be obtained 
because of the impossibility of applying the given boundary 
conditions along the odd-shaped boundary of a turbine 
blade (fig. 5) . A numerical solution is available, however, 
by application of the relaxation method (reference 5). A 
sketch of the blade cross section is covered by a network of 
points (a square network was chosen). Large net spacings 
are recommended at first, and the insertion of additional 
net points may be made at any time during the solution, 
thus permitting the use of a final network of any desired 
size. Temperature estimates for the net points can be obtained 
from the solutions of the chordwise blade-temperature- 
distribution equations for approximate shapes previously 
presented and from application of the boundary conditions 
given by equations (29) and (30).. Residuals, which may 
be, considered as interior heat sinks, can then be calculated 
at each net point whose immediately neighboring points 
remain witliin the boundary from the relation 

$o = 0i + 0s + 0a + 04 — 40o (31) 

where the subscript zero denotes a point in the blade cross 
section and subscripts 1, 2, 3, and 4 denote the points in 
the square surrounding the point with zero subscript. 



The object of the relaxation is to reduce the values of 
Q to zero, or as close to zero as possible. When the relaxa- 
tion equation (31) is employed^ the following procedure is 
used._ A change in 0 O alters Q u Q i; Q s , Q A by the same change 
and <2o by minus four times this change, all other values of 0 
remaining fixed. 

Equation (31) is the finite-difference equation correspond- 
ing to the partial differential equation (28). For net points, 
some of whose immediately neighboring points lie outside the 
boundary, equation (31) must be modified. For example, if 
point 1 lies outside the boundary and point 5 is the boundary 
point between 0 and 1, the following equation applies: 

$ 0 = 02 + 03 + 04 +^ — ( 32 ) 

where cf is the ratio of the distance between points 0 and 5 
to the net spacing. 

Corresponding changes in the values of Q result from the 
use of equation (32). No harm is done by overestimating 
the final values of 9, as successive calculations will establish 
them again. Continued relaxation eventually reduces all 
the residuals as desired and the blade temperatures can 
finally be obtained from the definition of 9, 

9=T t .,~T a 

APPLICATION OF ANALYSIS 

Typical numerical examples for a sample blade are 
presented to illustrate temperature trends and to determine 
the effects of various factors on the temperature distributions 
in a liquid-cooled turbine blade. A turbine blade whose 
external shape is similar to that of the rotor-root section of a 
conventional gas-turbine design was selected as the sample 
blade. Two 0.25-inch-diameter coolant passages were 
assumed in the blade (fig. 5), connected near the blade tip 
by a cross-over passage, and the cooling was assumed to 
occur tlirough forced convection. No allowance was made 
in these calculations for the effects of free-convection currents 
that might be present. A gas flow of 55 pounds per second 
(equivalent to a mass velocity of about 58 lb/(sec)(sq ft)) 
and a coolant flow of 7 pounds per minute per blade were 
assumed (to insure turbulent flow and to remove the depend- 
ence of tire blade-to-coolant heat-transfer coefficient on the 
length-to-diameter ratio of the coolant passage). 

For the one-dimensional spanwise case, in which the 
cooling was carried to the rotor hub, two 0.25-inch-diamcter 
coolant passages running circumferentially through the rim 
and ten 0.50-inch-diameter coolant passages running radially 
tlirough the rotor were assumed (fig. 1). For this case, 
water, ethylene glycol, and kerosene were all considered as 
possible coolants. Other calculations were made only for 
water as the coolant because of its superiority over the other 
coolants considered. 

Blade-to-coolant average heat-transfer coefficients were 
calculated by use of formulas presented in reference 6, 
page 168 for turbulent flow and page 190 for laminar flow 
(for ethylene glycol onl}'). Gas-to-blade average heat- 
transfer coefficients were calculated by use of formulas given 
on page- 236 of reference 6. The following coefficients were 
used in the analysis : 
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Outside heat-transfer coefficient, h, 

Btu / (hr) (sq ft) (°F) . _ 222 

Btu/(sec)(sq ft)(°F) 0.06167 

Inside heat-transfer coefficient, hi, t 
Water, 

Btu/(hr)(sqft)(°F) 2370 

Btu/ (seo) (sqft)(°F) 0. 6583 

Ethylene glycol, 

Btu/ (hr) (sqft)(°F) 649 

Btu/(sec)(sq ft)(°F) a 1803 

Kerosene, 

Btu/ (hr) (sq ft) (°F) 510 

Btu/(sec)(sq ft)(°F) 0.1417 

Other numerical values used in the calculations were 
( figs. I and 4): 

Effective gas temperature, T,,„ °F 2000-5000 

Average coolant temperature, Ti, °F 200 

Number of blades, Z 55 

Area of blade cross section 1, At, sqft 0. 00198 

Area of blade cross section 2, At, sq ft 0. 0013 

Area of cross section 3, Aj, sqft 0. 312 

Area of cross section 4, A*, sq ft _ — 0. 236 

Blade outside perimeter, l„ ft 0. 2542 

Blade inside perimeter, ft._ 0. 131 

Average radial distance, section 3, r Jt «„ ft. 0. 4917 

Average radial distance, section 4, r iia „ ft 0. 3333 

Length of chordwise trailing section, j, ft 0. 050 

Length of spanwise trailing section, 6, ft 0. 3333 . 

Thickness of trapezoid at coolant passage, r 2 , ft.. 0. 021 

Thickness of trapezoid at trailing edge, r t , ft 0. 003 

Thickness of rectangle, r, ft 0. 010 

Temperature at blade root, T r , °F 330 

Thermal conductivity, k B 

Btu/ (hr) (ft) (°F) 15-210 

Btu/ (sec) (ft) (°F)_ 0.00417-0.0583 

Cooling-air (in contact with rotor) temperature, T a , °F_ 0 

Heat-transfer coefficient between cooling air and rotor, 

Btu/(hr)(sq ft)(°F) 30 

Btu/ (sec) (sq ft)(°F) 0. 00874 

Average value in section 3 of h,ihij 
Water, 

Btu/ (hr) (ft) (°F) 30,360 

Btu/ (sec) (ft) (°F) 8. 433 

Ethylene glycol, 

Btu/ (hr) (ft) (°F) 8305 

Btu/ (sec) (ft) (°F) 2.307 

Kerosene, 

Btu/(hr)(ft) (°F) 6533 

Btu/(sec)(ft)(°F) 1.815 

Average value in section 4 of U,, 

Water, 

Btu/ (hr) (ft) (°F) 6107 

Btu/ (sec) (ft) (°F) L6E6 

Ethylene glycol, 

Btu/(hr)(ft)(°F) 1656 

Btu/ (sec) (ft) (°F) 0.460 

Kerosene, 

Btu/ (hr) (ft) (°F) 1318 

Btu/(sec)(ft)(°F) 0.366 


ONE-DIMENSIONAL SPANWISE TEMPERATURE DISTRIBUTION 

The constants of integration were found by solving equa- 
tions (5) to (10) for the specific examples stated; a blade 
thermal conductivity of 15 Btu/(br)(ft)(°F) was used. 
Blade coolant passages extending to within X« inch of the 
blade tip were considered in blades with spans ranging from 
lKe to 4Xs inches, with water as coolant, and for effective 
gas temperatures ranging from 2000° to 5000° F in order to 


study the effect of blade length on the temperature distri- 
bution. The t-emperature distributions were determined, 
by use of the calculated integration constants, from equa- 
tions (1) to (4) and are shown in figure 6. The general 
trend shown is a nearly constant temperature about equal 
to the coolant temperature through the rotor, a sharp tem- 
perature increase through the rim and the base of the blade, 
a nearly constant temperature (called the prevalent blade 
temperature) through the liquid-cooled part of the blade, 
and another sharp temperature increase to a value approach- 
ing the effective gas temperature in the uncooled part of the 
blade. From figure 6, it can readily be seen that the same 
prevalent blade temperature prevails through blades of 
various lengths cooled to within K« inch of the blade tip. 

In order to study the effect of coolant -passage length on 
the temperature distribution, other calculations were made, 
with water as coolant and for an effective gas temperature 
of 2000° F, for a 4Xs-mch blade span but with blade coolant 
passages of various lengths. The calculations, again deter- 
mined by use of equations (1) to (4), are plotted in figure 7; 
for short coolant passages high-temperature gradients exist 
throughout the blade, whereas for long coolant passages 
practically no temperature gradient exists throughout most 
of the blades. 



Fmcri 0.— Effect of various blade lengths on temperature distribution for water-cooled 
gas turbine for effective gas temperatures from 2000° to 5000° F. One-dimensional spanwise 
analysis. Coolant passages extend to within H« Inch of blade tip; coolant flow, 7 pounds 
per minute per blade; thermal conductivity of blade, 15 Btu/ (hr) (ft) (”F). 



Figure 7.— Effect of various coolant-passage lengths on temperature distribution. One- 
dlmenslonal spanwise analysis. Blade length, 4Jde inches; coolant flow, 7 pounds per 
minute per blade; thermal conductivity of blade, 15 Btu/(hr) (ft) (°F). 
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Calculations were also made for the same blades with 
4-inch coolant passages and for an effective gas temperature 
of 2000° F with water, ethylene glycol, and kerosene con- 
sidered as possible coolants. These results showing the 
effect of various coolants on temperature distribution are 
given in figure 8. For an effective gas temperature of 2000° F 
and a coolant flow of 7 pounds per minute per blade, the 
prevalent blade temperature in degrees Fahrenheit is about 
one fourth, two fifths, and one half of the effective gas 
temperature for water, ethylene glycol, and kerosene cool- 
ants, respectively. The consideration of cooling air on the 
outside of the rotor accounts for the rotor temperature being 
less than the coolant temperature. 

Finally, for the spanwise case, figure 9 shows the variation 
of coolant flow on prevalent blade temperature. The 
prevalent blade temperature decreases with increasing 
coolant flow ; the rate of this decrease diminishes as the cool- 
ant flow increases. It can be seen that the cooling effective- 
ness, defined as changes from 0.68 to 

0.91 for water as the coolant flow is increased from 2 to 16 
pounds per minute per blade. The superiority of water is 
apparent. 
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Figure 8.— Temperature distribution for gas turbine for various liquid coolants. One- 
dimensional spanwise analysis. Coolant-passage length, 4 Inches; coolant passages extend 
to within )i« inch of blade tip; coolant flow, 7 pounds per minute per blade; thermal con- 
ductivity of blade, 15 Bto/(hr) (ft) (°F). 



Figure 9.— Variation of prevalent blade temperature with coolant flow. One-dimensional 
spanwise analysts. Thermal conductivity of blade, 15 Btu/(hr) (ft) (°F). 


THREE-DIMENSIONAL TEMPERATURE DISTRIBUTION THROUGH 
— SIMPLIFIED TRAILING SECTION 

A three-dimensional temperature distribution in a turbine- 
blade trailing section approximated by a rectangular paral- 
lelepiped (fig. 2) was determined by use of equation (12) 
(approximated by use of m=l only, as previously mentioned). 
The constants L n , Q ni and TV) were obtained from equations 
(16), (17), and (18), respectively; the constants K Un and 
0 Un from equations (14) and (15), respectively; and the 
constants M lt „ and P i,„ from equations (13). The distribu- 
tion was found in two planes representing the maximum and 
minimum temperatures for the Z-axis ; the first plane was 
located, at the side of the cross section of the rectangular 
parallelepiped and the second plane was on the median plane 
through the rectangular parallelepiped. The results are 
shown in figure 10(a) for an effective gas temperature of 
2000° F and a thermal conductivity of 15 Btu/(hr)(ft)(°F) 
and give the temperature distribution at various distances 
from the coolant passage. The curve labeled “approxima- 
tion of temperature” in figure 10(a) is a one-dimensional 
chordwise distribution through the approximated rectan- 
gular trailing section. Similar results, for thermal conduc- 
tivities of 120 and 210 Btu/(br)(ft)(°F), are given in figures 
10(b) and 10(c), respective^'. Calculation of cosTVir/2 (see 
equation (12)) reveals the temperature variation in the. two 
planes to be about 3.6 percent of 8 for a blade thermal con- 
ductivity of 15 Btu/(hr) (ft)(°F), 0.5 percent for a thermal 
conductivity of 120 Btu/ (hr) (ft) (°F), and 0.3 percent for "a 
thermal conductivity of 210 Btu / (hr) (ft) (°F) . For a thermal 
conductivity of 15 Btu/ (hr) (ft) (°F), figure 10(a) shows a 
constant temperature for the last three quarters of the blade; 



(a) Thermal conductivity of blade, 15 Btu/(hr) (ft) (°F). 

Figure 10.— Three^imensional temperature distribution in trailing section of turbine blade. 
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Blade length, in. 

(b) Thermal conductivity of blade, 120 Btn/(hr) (ft) (°F). 

FlilfRS 10.— Continued. Three-dimensional temperature distribution In trailing section 

of turbine blade. 



Fi'iURE 10— Concluded. Three-dimensional temperature distribution In trailing section 

of turbine blade. 



Fiscal 11. — Comparison of one -dimensional chordwlse temperature distribution through 
rectangular trailing section with three-dimensional temperature distribution through rec- 
tangular parallelepiped. Thermal conductivity of blade, IS Btu/(hr) (ft) (°F). 


that is, conduction to the rim affects about the first quarter 
of the blade length. Figure 10(a) also shows that the level 
of the temperature in the principal portion of the blade rises 
rapidly as the distance from the coolant passage is increased. 
Near the coolant passage, the distribution is in good agree- 
ment with the one-dimensional spanwise distribution pre- 
sented in figures 6 and 8. As the thermal conductivity is 
increased, as shown in figures 10(b) and 10(c), the part of 
the blade affected by rim conduction increases; for a thermal 
conductivity of 210 Btu/(hr)(ft)(°F), about two thirds of the 
blade length shows this effect. 

ONE-DIMENSIONAL CHORD WISE TEMPERATURE DISTRIBUTION 
THROUGH SIMPLIFIED SHAPES 

Because the three-dimensional distribution resulted in a 
constant spanwise blade temperature in the region of the 
blade bej’ond the influence of rim cooling and because the 
critical blade point, as determined from a temperature-stress 
relation, may likewise be beyond the region of run-cooling 
influence, one-dimensional chordwise temperature distribu- 
tions were obtained. In order to compare the three-dimen- 
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Fujurb 12— One-dimensional chord wise analysis of effect of shape on prevalent blade tem- 
perature. Thermal conductivity of blade, 15 Btu/(hr) (ft) (°F). 

sional distribution with, a one-dimensional chordwise 
distribution, a rectangular trailing section was first consid- 
ered. Figure 11 shows this comparison for an effective gas 
temperature of 2000° F and a blade thermal conductivity of 
15 Btu/(hr) (ft) (°F) ; the one-dimensional distribution was 
determined by use of equation (19). From figure 11 it can 
be seen that the trailing-edge temperature increases from 
1500° to 1850° F as the distance from the blade root is 
increased from 0.25 to 1.0 inch for the tlu'ee-dimensional case, 
even for low-conductivity material. On the other hand, in 
the region of the blade beyond the influence of rim cooling 
(when the distance from the blade root is 1 in. or more), the 
one-dimensional and three-dimensional results approach 
identity; for such a region, a three-dimensional solution is 
unnecessary when low-conductivity materials are considered. 

In order to more nearly approximate a blade trailing sec- 
tion, a trapezoidal approximation was considered. For 
comparative purposes the trapezoidal and rectangular sec- 


tions were constructed to have equal lengths and areas. In 
general, however, the trapezoid is constructed so that the 
thickness at the coolant passages equals one-half the actual 
blade surface exposed to the coolant. This dimension, in 
turn, fixes the length of the trapezoid. 

The temperature distribution through a trapezoidal sec- 
tion was obtained by use of equation (20). Such a distri- 
bution, for an effective gas temperature of 2000° F and a 
blade thermal conductivity of 15 Btu/(hr)(ft)(°F), is 
compared with the distribution tlrrougli a rectangular sec- 
tion in figure 12. The temperature distribution for the 
trapezoidal section has a slightly steeper slope than that for 
the rectangular section at distances remote from the coolant, 
passagesyand the temperatures at the trailing edge and the 
coolant-passage wall are lower. Part of this lower tempera- 
ture forthe trapezoidal section is due to the additional thick- 
ness of tire cooling surface. 

The effect of varying thermal conductivity from 15 to 210 
Btu/(hr).(ft)(°F) on the temperature distributions, along the 
center line of rectangular and trapezoidal sections is shown 
in figures 13(a) and 13(b), respectively. In each case, as 
thermal conductivity is increased, the temperature distri- 
butions flatten and approach linearity, decreasing in the 



(a) Rectangular trailing section. 

Figure 13.— rOne dimensional chordwise analysis of prevalent blade temperature for trailing 

section of blade. 








w Coolant temperature 

200 ' 1 1 1 — < 

0 J .2 .3 A .5 .6 

Distance from coolant passage, in. 

(bi Trapezoidal trailing section. 

Fioi’RB 13.— Concluded. One-dimensional chordwtse analysis of prevalent blade tempera- 
ture for traHln* section of blade. 

trailing section and increasing near the coolant passages. 
Figure 14 shows blade trailing-edge temperatures for various 
thermal conductivities obtained by use of the rectangular 
and trapezoidal approximations. Cooling is substantially 
greater for the more representative trapezoidal section. For 
currently used high-temperature alloys with thermal 
conductivities in the neighborhood of 15 Btu/(hr)(ft)(°F), 
figure 14 shows a trailing-edge temperature difference for 
the rectangular and trapezoidal sections of only about 60°F 
for an effective gas temperature of 2000° F. 

Chordwise temperatures were also calculated, by use of 
equations (22) to (27), for concentric circle annuli and 
parallel-plate approximations. These temperatures are not 
plotted in this report; they were only used as initial approxi- 
mations for the two-dimensional numerical calculations made 
for an actual blade shape. 

TWO-DIMENSIONAL CHOBDWISE TEMPERATURE DISTRIBUTION 

Two-dimensional temperature distributions were deter- 
mined for an actual blade shape by application of the relaxa- 
tion method. Separate calculations were made for the cases 
where an average gas-to-blade heat-transfer coefficient and 
where a typical variation in the gas-to-blade heat-transfer 
coefficient as shown in figure 15 were considered. Initial 


Gos-fo-b/ode heof- 
tronsfer coefficient 
(Bfu/(hr) (sq ft) (°F)) 
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Figckz 15 —Cross section of turbine blade showing typical variation of beat-transfer coeffi 
dent around blade perimeter. 
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trial solutions were determined by use of equation (20) for 
the trailing section, equations (22) to (24) for the leading 
section, and equations (26) and (27) for the thin-wall sections 
of the blade near the coolant passages. A comparison of 
the temperature distributions for assumed variable and 
constant gas-to-blade heat-transfer coefficients, for an 
effective gas temperature of 2000° F, an average coolant 
temperature of 200° F, and a blade thermal conductivity of 
15 Btu/(hr)(ft)(°F) is shown in figure 16. The blade 
temperatures obtained are nearly the same for both cases 
except at the leading and trailing edges. Use of the average 
coefficient gives a conservative estimate of the trailing-edge 
temperature and a temperature that is somewhat too low 
near the leading edge. 

A calculation has also been made for the average gas-to- 
blade coefficient and for a blade with a thermal conductivity of 
100 Btu/(hr) (ft) (°F). A comparison of this solution with the 
similar one for a thermal conductivity of 15 Btu/(hr)(ft) (°F) 
is shown in figure 17. The high-conductivity blade has 


Gas-to-b/ode 

heat-transfer 

coefficient 



Figure 16.— Effect of variation of gas-to-blade heat-transfer coefficient on temperature 
distribution through crosa section of water-cooled turbine blade. Thermal conductivity 
of blade, 15 Btu/(hr) (ft) (°F); effective gas temperature, 2000° F; average water tempera- 
ture, 200° F. 


Thermo ! 
conductivity 



Figure 17. — Effect of variation of thermal conductivity on temperature distribution through 
water-cooled turbine blade. Gas-to-blade heat-transfer coefficient, 222 Btu/(hr)(sq ftK'Fj; 
effective gas temperature ,2000° F -average water temperature, 200° F. 


about a 250° and a 600° F lower temperature than the low- 
conductivity blade at the leading and trailing edges, respec- 
tively. Little temperature difference is obtained in the 
center of the blade, where extremely good cooling prevails. 

The results of the various relaxation solutions are in good 
agreement with the one-dimensional chordwise approxima- 
tions that were used to start the numerical solutions; that 
is, good, representative temperatures are obtainable by use 
of one-dimensional chordwise approximations. 

In order to determine the effect of distance from the cool- 
ant passage on the trailing-edge temperature, trailing-edge 
temperatures were determined for the liquid-cooled blade 
with five coolant passages shown in figure 18. Various 
length trailing sections were obtained by successively reduc- 
ing the length of the trailing section shown in figure 18. The 
temperatures were determined by use of equation (20) and 
are shown for various thermal conductivities in figure 19. 
The trailing-edge temperature is reduced almost linearly as 
the length of the trailing section is decreased. The effect of 
thermal conductivity also decreases as the trading-section 
length is decreased. 



Figure IS.— Cross section of water-cooled turbine blade showing location and slr.e of live 
coolant passages UBed In determination of effect of varying tralllng-scctlon length on 
trailing-edge temperature. 


NONDIMENSIONAL CHARTS 


The availability of several nondimensional charts, to be 
subsequently discussed, eliminates the necessity for some 
numerical calculation. 

The prevalent blade temperature is given by equation 
(2a) as 


B - r ~ hi'+KA* 


■ (2a) 


After division by A-,, 2 / 1 , 2 , this equation may be written 


where 


rp Tj-\-\T s , t 

B ' v 1+X 


X 


At, zh. s 
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Fiocri 19.— TratHng-edge temperature as function of length of trailing section of blade 
for water-cooled turbine blades with five coolant passages and various thermal conduc- 
tivities. 


Subtraction of T f _ t from both members of this equation 
leads to 

Tg, t — T B ,p 1 

1+X 

y y 

The plot of against X is shown in figure 20. 

1 1, e — i 1 

For any given blade, X can be evaluated and the value of 

rp p 

~ - can be obtained. Finally, for the desired effective 

l(,i -t i 

gas and coolant temperatures, a single simple algebraic oper- 
ation results in the desired value of the prevalent blade 
temperature. 

Another nondimensional chart, which gives the one- 
dimensional chordwise temperature distribution through a 
rectangular section, is also presented in figure 21. The 
temperature distribution through a rectangular section is 
given by equation (19), 

^ (TV.- TO cosh <py' 

— —j (i9) 

<p sinh <pj ,J rr~ cosh <pj' 

h*B 


1S1 



Fiocrx 20. — Nondimensional chart for use in determination of prevalent blade temperature. 


For a given turbine and set of turbine operating conditions, 
values can be determined for all the quantities in this equa- 
tion. A semilog plot with <pj r as abscissa and ht/(pk B as a 
parameter results in a nondimensional cooling ratio 


y y s 

(Tg.g—Ti) cosh <py' 


The addition of a second quadrant, with cosh <py' as param- 

y y 

eter, yields values of the temperature ratio — Tfr at any 

point in the rectangular section. This chart is given in 
figure 21. The chart is used as follows : A vertical line is 
constructed through the calculated value of the abscissa vj' 
and is extended to the calculated value of the parameter hil<pk B . 
From this point, a horizontal line is drawn, extending into 
the second quadrant, and intersecting several lines repre- 
senting various values of cosh <py' . Vertical lines from the 
intersection points to the abscissa in the second quadrant 

y y 

give values of the temperature ratio ™*- — sft at various 

■i t . « i 


positions in the rectangular section. Values of the tem- 
perature T b are then easily obtainable. 

In order to illustrate the use of figure 21, a rectangular 
section with dimensions equal to those previously tabulated 
will be considered. For a blade thermal conductivity of 
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Floi'KX 21.— Nondlmenslona] cliart for use In determination of temperature distribution through rectangular trailing section of cooled turbine blade. 


15 Btu/(hr)(ft)(°F) and outside and inside heat-transfer 
coefficients of 222 and 2370 Btu/(hr)(sq ft)(°F), respec- 
tively, the parameters used in figure 21 are evaluated as 
follows: 

j h 2370 

<f>k B (54.4) (15) a 

The point represented by these parameters is given by the 
circle in the first quadrant of figure 21; the square on fig- 
ure 21 designates the corresponding point for a thermal con- 
ductivity of 210 Btu/(hr) (ft)(°F). A horizontal line pass- 
ing through the circle and intersecting the family of lines in 
the second quadrant of figure 21 gives as the abscissa of the 

y y 

second quadrant the values of the temperature ratio — sr 

4 1 , e — 4 l 

for various positions in the rectangular section as follows : 

£*=0.075 for W '= 0 
= 0.115 for <py' = 1 
= 0.285 for <py'=2 
= 0.75 for <py'=3 

For 2Y»=2000° F and 7 1 i =200° F, it follows that the blade- 
temperature calculation results in the following values: 

T fl =1865 for <py'=0 (trailing-edge temperature) 

= 1793 for <py'=l 
= 1487 for ^'=2 

=650 for <py'— 3 (temperature at coolant passage) 


These results compare favorably with the calculated distri- 
bution shown on figure 1 2. 

CONCLUDING REMARKS 

Analyses have been presented for obtaining spanwise 
temperature distributions near a coolant passage, chordwise 
distributions (for both approximated and actual shapes) in 
regions where conduction to the run is inappreciable, and 
three-dimensional distributions for approximated blade 
shapes. Numerical examples based on specific blade con- 
figurations and heat-transfer coefficients available from 
unclassified sources have been presented. Although the 
analyses are exact, the numerical values used in tire calcu- 
lations may not necessarily be the same as those for cooled 
turbines. The numerical examples have been presented to 
indicate the range of applicability of the various analyses 
and to present the general nature of temperature distribu- 
tions in liquid-cooled turbine blades. In the following 
paragraphs the more important results are reviewed to show 
their general guidance for design studies. 

The three-dimensional temperature distribution includes 
conduction to the rim, whereas the two-dimensional and 
one-dimensional chordwise distributions do not. Conse- 
quently, the temperatures obtained from a three-dimensional 
investigation are less than those otherwise obtained. As a 
result, it is advisable to determine the simplified solutions 
first; if the temperatures that result are not excessively high, 
a three-dimensional investigation is unnecessary; whereas if 
the resulting temperatures are large, a three-dimensional 
investigation can be made. 
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The two-dimensional solution just referred to is of necessity 
a numerical solution (because of the varying boundary con- 
dition caused by the shape of the blade). It has been shown 
by the calculations previously presented that for uniform 
outside heat-transfer coefficients, simplified one-dimensional 
distributions were in excellent, agreement with the two- 
dimensional relaxation solution. When variable outside 
heat-transfer coefficients were considered, the simplified 
solutions indicated optimistic results and the two- 
dimensional relaxation solution appears essential. Increasing 
blade thermal conductivity results in raising the coo ling 
surface temperature and lowering the trailing-edge tempera- 
ture and thus in a more uniform blade temperature. 


One-dimensional spanwise distributions proved to be valid 
near the coolant passages. The temperature of the cooled 
part of the blade (prevalent blade temperature) was inde- 
pendent of blade length. If the inside heat-transfer coeffi- 
cient is increased, the prevalent blade temperature decreases; 
it was shown for a particular case that the prevalent blade 
temperature for water as coolant was only about half that 
for kerosene as coolant. 

Lewis Flight Propulsion Laboratory 
National Advisory Committee for Aeronautics 
Cleveland, Ohio, October 27, 1950 


APPENDIX A 


The following symbols are used in this report: 
area of metal, sq ft 
blade span, ft 

., Cu constants of integration 

diameter of circle used to approximate leading 
section of blade, ft 

ratio of distance between net points 0 and 5 
and net spacing 

Hankel functions (special kinds of Bessel 
function) 

heat-transfer coefficient, Btu/(see)(sq ft)(°F) 
or Btu/(hr)(sq ft)(°F) 

Bessel functions 

chordwise distance from blade trailing edge to 
coolant passage, ft 

( — V 

V k B sin 

thermal conductivity, Btu/(sec)(ft)(°F) or 
Btu/(hr)(ft)(°F) 
perimeter, ft 

0, P points on figure 3 
heat flow, Btu/sec 

radial distance from hub of turbine, ft 
temperature, °F 

distance from blade tip to blade element, ft 
distance from trailing edge to blade element, ft 
number of blades 

distance from median plane of section to blade 
element, ft 

(t f il^oftfi) i i (if i )HI» (if a) 


.1 
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<u 
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SYMBOLS 

Subscripts: 
a 
av 
B 
e 

9 
i 
l 


2K 


~ , . r. q— tanjf B 
_ ■ 2 tan \p 


bo Ip 
ht a 1 j, 2 

distance between parallel plates, ft 
thickness of trailing section, ft 

(fey 


Superscript-: 

prime 

Functions: 


B 


D 


E 


-"feif) 


air 

average 

blade 

effective 

gas 

inside (inner with D) 
liquid 

summation indices 

outside (outer with D) 

prevalent 

rotor 

blade root 
rim 

blade tip 

denote sections in spanwise investigation 
when used with*; 1 and 2 denote trapezoidal 
thicknesses at trailing edge and coolant 
passage, respectively, when used with r 
denote end near blade tip and end near rotor 
hub of various sections 

linear dimension increased by r/2 

( KU Y 

\k B A Bt iJ 
f Klo-kht'jli't S. i 

\ k B A b , 2 J 

\ kgAs.t ) 
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K, L, M, N, integration constants (See equations (13) to 
0, P, Q (18).) 

R k B M m , n sink M m , n j'+h t cosh M n j' 

4k . f — 11“— 1 


2(T t , e -T B ,r ) 



sin 2 Q n y 

~2Qut: 


The heat balance for the uncooled, section of the blade 
(fig. 1, section 1) is given as follows: 

d IJ 1 

Heat entering by radial conduction= —k B A B , i 1 

Heat entering sides by convection ^^(TY,,— T BA )dxx 
Heat leaving by radial conduction 

=-j kaA *- 1 Jr, ( Tb ’ dXi ) 

The heat-balance equation is 

- k B A B , , + K lo(T,. e~T BA )dXi 

=~ k * A *^( TB -'+ , ^ dx ') 


APPENDIX B 

ONE-DIMENSIONAL SPANWISE TEMPERATURE DISTRIBUTION 

where C 3 and C\ are integration constants and 


or 


d*T B , , 

dxS 


where 

A solution is 


■A i T B ,x=—A i T t>[ 


I / M» V 
\k g A B . 1/ 


(Bl) 


T BA =T,',— Ci cosh A(zi+ Q 


( 1 ) 


where C\ and C s are integration constants. 

For the cooled section of the blade (fig. 1 , section 2) . the 
heat balance is given as follows: 

dT 

Heat entering by radial conduction =—k B A BA 
Heat entering sides by convection=AX(?Y,— T B ^dx t 
Heat leaving by radial conduction = 


—k B A 


B*±B,2 


dx 


dx 2 ) 


Heat leaving by radial cony&ction=h i- J,, A {T BA —Ti) dx 2 
The heat-balance equation is 
dT B 


-k B A 


B, 3 ‘ 


dxt 


^+hJ 0 (T^,-T B , 1 )dx 2 


— ~^k B A B ' j dXi^-\-ht A lt A {T Bt i — Ti) dx a 


or 

where 

and 

A solution is 


d*T, 

dx« 2 


b ' 2 -B 2 T b , 2 =—D- 


(B2) 


B 


/ hi ~|~ hj t a lj t 

\ k B A B2 ) 


D=(± 


h„h Tt, e ~r hi' 1 


k B A Bi a ) 

T BA =T B , v -C % e B *t- C&~ B h 


( 2 ) 


it haloT fi ,-\-ht'J,i'iTi 

Xb -’-\b)- Kl 0 +h t ji, t 

For the rim section (assumed constant area) of the rotor 
(fig. 1, section 3), curvature was neglected because the rim 
thickness was small in comparison with the rim radius. In 
addition, average values of h. iA l (A and r s were used. The 
heat-balance equation reduced to 


d s T r i 


E‘T^=-F* 


(B3) 


where 


and 


E--- 


_ 4Tr 8| at h a ~f~ (ht, di, s)a 

k rlmA r i m 


p2 4»T s, as kg Pg ~[~ (Af ; di, s)agTf 

k t i m A tim 


A solution is 


(31 


T T i m =G+ C t e Bl * + C 6 e _Sl » 

where C s and C* are integration constants and where 


For the rotor section (assumed constant strength) (fig. 1, 
section 4), average values are used for and for r 4 , and 

an approximate solution (as in section 3) is obtained. Be- 
cause the rotor is liquid-cooled, such a solution is adequate. 
The resulting differential equation for this section is 


d 2 T R 


-WT^-F 


where 


and 


P- 


dr* 

rj 4xr 4 t aihgA" (Af, 4 Y 4 )flp 

k R A R 

^■xr^a,haT a -{-{hi A li^ a ,Ti 

k R A R 


(B4) 


A solution is 


T r =J + C) cosh H r 

where C 2 is an integration constant and where 

/V 


( 4 ) 


The boundary condition 

dTg 

dr 

has already been applied. 


HU 


= 0 for r=0 
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A three-dimensional temperature distribution through a 
rectangular parallelepiped, an approximation for the trailing 
section of a liquid-cooled turbine blade, is given as follows 
(fig. 2): 

dT B 


Heat entering element from to^= — k B dydz 


da: 


dT B 

by 


Heat entering element from right end = —k B dxdz 

dT 

Heat entering element from front =—k B dxdy 

Heat leaving element at bottom 

, . . /bT B , d 2 T s , \ 

= - w ' <I (a»- + w iz ) 

Heat leaving element at left end 

--*.*»* (^+^*i») 

Heat leaving element at rear 


=—k B dxdy ( 


i>T B d a r B 
, dz + dz 2 


d z 


) 


The heat-balance equation is 


i j j &T B i j j d T b r 7j dlZa 
— k B dydz —k B dxdz k B dxdy ~^r= 


— kadydz (^+*?*dx)-k B dxdz dy)- 

k B dxdy ^+^ s d 2 ) 

d *T B , d 2 T B , d 2 T fl n 


dr 2 dy* d z 1 


(Cl) 


Simplification in applying boundary conditions (as fully 
described on p. 4) and use of the effective gas temperature 
T t , t instead of the blade temperature T B as the reference 
temperature, that is, use of the substitution 

d=T t ,.-T B 

< 

change equation (Cl) into 


d ! 0 d 2 0 d ! 0 

bx'-dy'^dz*— 


( 11 ) 


The origin of the coordinates chosen is shown hi figure 2; 
the plane 2=0 is the median plane of the rectangular para- 
lelepiped and from considerations of symmetry no heat flows 
across this median plane. The boundary conditions to be 
applied are 


S^,=0 when r'=0 
ox' 

(C2) 

d0 

fy,= 0 when y'=0 

(C3) 

5— =0 when 2 = 0 
oz 

(04) 

d=T t e — T b>t when x'=b' 

(C5) 


and 


ks Ey =ht {T t ,,—Ti—d) when y'=j' (C6) 


ks^=— A„0 when z=r/2 


(C7) 


A solution of equation (11) satisfying the three boundary 
conditions (C2), (C3), and (C4) is 


CD 00 

S=DS (K» t , cos L„x' cosh M n , H y' cos N M z + 

m— 1 «— 1 

O m , n cosh P m ,*x' cos Q x y' cos N m z) (12) 

'where K, L, M, N, 0, P, and Q are constants. Relations 
among these constants are 


M m .S=N m *+LS) 

*=#»*+ «. 9 J 


(m,n= 1,2,3, 




(13) 


The boundary^ condition expressed by equation (C5) leads 
to a determination of L. Substitution of equation (12) into 
equation (C5) gives 


T,.,— T B , r = S X) i K *.* cos L n b' cosh M m , % y' cos N n z-\- 

m— I n-1 

0 m , » cosh P„, n b' cos Q k y' cos N m z) (C8) 

It can be seen that the first member in. the double s umma tion 
in equation (12) vanishes when cos L A x' — 0; hence, if 



from equation (C'8) it follows that- 


Tt,i—T B , r = X X °m,m cosh P m , , b' cos Q H y' cos NnZ (C9) 

n»— I n— I 

This is a Fourier development along y' and z. 

The development of equation (C6) leads to a determination 
of Q n . Equation (C6) becomes, with the aid of equation (12) 
and its derivative, 


X X £b[A„,„(cos L k x 0 M„, a sinh cos N m z— 

m-I n-1 

O m , n (cosh Pn r nX')Q n sin Q a j' COS N nZ\ 

= X X htiTte—Tt—Ka^cvs L n x' cosh M u ,J'coaN m z- 

m—l n-1 

O n ,n cosh Pn.nX' COS Q„ / COS N m z) (ClO) 

In order to simplify equation (ClO) at the boundary when 
y'=j' and to solve for the constant Q m , the second summation 
in the left member is equated to the last summation in the 
right member term by term. Then 

^ <(T t , — r,) = x x A 'n. «(k B M m , „ sinh M n , J' + 

ra-1 n-1 

ht cosh M n , n f) cos L*x' cos N m z (Cll) 
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2 cosh P n , a x')Q n sin Q n j' cos N n z 

m- 1 n-1 

= 2 cosll P m , n x' COS Q n j' cos N m z (C12) 

1 n— 1 

The sums are equated term by term. Therefore, Q n is 
determined as any solution of ... 

hif_ 

hi kg 


tan (QJ')= 




(17) 


The development of equation (C7) leads to a determina- 
tion of N m . 

2 2— cos LnX' cosh M„,„y' (lV m sin N n £)+ 
0„ t n cosh P m , a x' cos Q n y' (N m sin N n £Y|=2 2“ 

\ 4/J 771-1 71-1 

h„ (Km. n COS LnX' COsh M n , nV' COS N„ ~ + 

• 0 Mi n cosh P m _ nX' cos Q n y' cos N m 0 (Cl 3) 

or 

h 0 r 


. T Kq 2&5 

t&nNm 2-kjr m ~rn^ 


(18) 


The possibility of determining values of K mt , and O m% K to 
satisfy equations (C9) and (Cll) has been established in 
textbooks on Fourier and other harmonic series (reference 7, 
pp. 118-121), and it is only necessary that the values 
determined define a convergent series. 

Values of K n „ are determined by integrating equation (Cll) 
between the limits x'=0 to x'=b' and 2=0 to z=r/2 and 
substituting the values previously determined for L n and 
N„ in equations (16) and (18). The integration is accom- 
plished in two steps using the functions cos L,x' dx' and 
cos N,z dz as multipliers, where $ and v are integers. For 
abbreviation, 


Then 


R=k B M m , n sinli M mtn j'+h t cosh M n , n j' (C14) 


rx'-b’ 

» cos N m z cos 2 L n z r dx' 

Jz'-O 


m-1 


because all the terms 


=h i (T 3 ' t —T l ) f cos L n x' dx' 

Jx'-O 
ms 

I cos L„x' cos L,x' dx' 

Jo 


(Cl 5) 


or 




cos .V ,2 


sin Ln h’ 


—2h t (T lte —Ti) 


When equation (16) is used for L„, there results 

2 RK mi n cos N m z=ih t (T t .'-T l ) 

m - 1 * [Xn—i) 

When equation (Cl8) is integrated in terms of z, 
the multiplier cos M, z ds, it follows that 


f‘~2 /**"2 

RK„ t1i cos 2 N m z dz=Sn cos N n zdz 
t ’ Jt-0 J«- 0 

because all the terms 

r 

J* cos N m z cos N,z dz . 
vanish if m^v. Upon integration, 

RK m n / . sin 2V m z V S n / . Ar \* ' 


2JV. 


or 


RK, 


l sin 2N m |\ 2S n sin N m 
l ' n \i + 2 N„ / NnT 


Therefore, 


K„ 


2 Sn 8iQiVw i 1 

sin 2N m ^ 

\ 2iV " \ J 


Similarly, from equation (C9) it is found that 


sin N m 


O m .n=2T n 




cosh I 1 


Lu^i) 

\ 2i¥ -§ j 


where 


nn 9 (rp rp \ SID. Qnj t 

l B , r ) QJ , s [ n2 Q u j' \ 

2 Qnj' ) 


(Cl 6) 

(Cl 7) 

(CIS) 

using 

(C19) 

(C20) 

( 021 ). 

(C22) 

(C23) 

(C24) 


vanish if s^n. Integration of equation (C15) and sub- 
stitution of limits lead to the result 
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Finally, substitution of equations (C14) and (Cl8) in equation (C22) and substitution of equation (C24) in equation (C23) 
lead to the following values of- the coefficients K m * and 0 m , *: 


and 
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Rectangular trailing section of blade removed from in- 
fluence of rim cooling. — The heat balance for an element of 
the rectangular section is given as follows (fig. 4(a)): 

Heat entering by conduction (right end )=k B r-^ 

Heat entering by convection (sides) =2 h 0 9dy 
Heat leaving by conduction (left end) 

The heat-balance equation is 

k.r'^+2hMy=kn,(^+ d ^ i v) 


or 


where 


A solution is 




( 2K\i 

^-yjTr) 

8=C S cosh <p{y+C 9 ) 


(Dl) 


(D2) 


where C« and C\ are integration constants. The boundary 
conditions to be applied in the evaluation of the integration 
constants are 


k BT , lr=h i T{TB-T l )=hMT t '.-T l -9) when y=j (D3) 

ay 


and 


k B ^=h 0 d when y=0 


(D4) 


From equations (D2) and (D4), it follows that 
k B C\<p sinh (<pCg)=h 0 C s cosh (#>C e ) 


a 




tanli" 




<p 


(D5) 


when r is not too large. Calculations revealed that the. 

difference between C a and r/2 varied from 0.4 to 0.5 percent 
for the values of <p considered and had no appreciable effect 
0.005 foot or more inside the trailing edge. 

From equations (D2), (D3), and (D5) it follows that 
k B rC s <p sinh (<pj')=h i T[T t ,.—T,—C» cosh 
and therefore 






'S — h 

• <p sinh (<pj')+y cos h (‘Pj') 


(D.6) 


Substitution of equations (D5) and (D6) in equation (D2) 
leads to the final equation 


0= 


jr (T^—Ti) cosh (ipy') 

<p sinh cosh (<pj r ) 


(19) 


Trapezoidal trailing section of a blade removed from in- 
fluence of rim cooling. — In reference 8, the temperature- 
distribution equation along the axis of a' trapezoidal section 
is derived. The solution obtained is expressed by the rela- 
tion 

e=C l0 J o (it:)+CniH o (it) (D7) 

where (?i 0 and C n are integration constants, J„ and iH„ are 
Bessel functions, and f is defined as 


r*=4K 2 y' 


+ 


n(l — tan $) 
2 tan 


] 


(D8) 


and therefore 
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where 

and 
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K 2 -- 


ka 


k B sin 


(D9) 


(DIO) 


The constants C ir> and C n in equation (D7) may be evaluat ed 
by application of the following boundary conditions: 


and 


or 


and 


^=0 when y '= 0 


do 


k B ~~=h t (T t , t -T t -B) when y’=j' 
when Jr=fr 


k B 2 K 1 
de- 


'■kiiT^t—Ti—S) when {-=f 2 


It follows from equation (D8) that when y'= 0, 


1 2/ L 2 tan $ J 


and when y'—j', 


t- i r i(l tan ^)“| 1 

f ’- 2 K[J + 2 tan $ J 


2 tan ^ 

Differentiation of equation (D8) gives 


hence, 


df 2k*dy' 

de n . de n 

t-= 0 when -j— ,=0 

dr dy' 


From the properties of the Bessel functions, 
dJ„(y)_ 


and 


dy 

dH 0 (y) 


■ Ji(y) 


dy 




Differentation of Equation (D7) therefore gives 

cwji(tr)+c'„n i (ift 

and it follows that 


Cn- 


: ZTiCtfi) 


d0 

dr 


=0 


(Dll) 
(D12) 
(D 13) 
(D14) 


(D15) 


(D16) 


(D17) 


(D18) 


(D19) 


From equation (D12), with the use of equations (D7) and 
(D17), it is found that 


*i 


2 ff 2 

r 2 


t-c\ 0 t«/ 1 (ir J )+c , 11 //,(ir 2 )] 


— — Ti — r 2 ) — C n iH 0 (i r 2 )] (D20) 


The values of the integration constants Cm and Cn are 
now found by solving simultaneously equations (D19) and 
(D20). Insertion of these values in equation (D7) gives the 
desired equation for the temperature distribution through 
a trapezoidal section 


where 


2^ iT s , e -T t ) [H,(»ri) ^(irt+i^iCtrowtr)] 
U J x a rofli« r*)] - [ffi(i rot «*(* r*)] + 2 ^f- r 


r=ff 1 (iri)J»(irs!)+id'i(iri)td4(ir*) 


( 20 ) 


For a wedge, r^O. Then, from equation (D 15) it is seen 
that r i — 0, and as a consequence, //lCifi)* ro and it/i(ffi)=0. 
Equation (20) then simplifies to 


2 ^ ^ro-^.(iro 


( 21 ) 


Concentric circle annulus approximation for blade leading 
edge. — The derivations of the equations giving blade tem- 
peratures at the inner and outer edges of the annular region 
are given as follows (fig. 4(c)): 

Heat entering from hot gas=2ir^ h 0 (T,, t — T B ,o) (D21) 


Heat flowing over circle with diameter 
D= 2t ®k B - dT * 


(!) 

D t 


(D22) 


Heat . leaving to coolant=2ir h t (T B , t — T t ) (D23) 


A solution of equation (D22) is 

t -=(Yk)'A c 4) 


(D24) 


where C x2 is an integration constant. Application of the 
boundary conditions 


and 


Tb=T b , 0 when D = D„ 
T B =T Bri when D = D t 


to equation (D24) results in the elimination of C xi and an 
evaluation of Q, namely, 

r D 0 




when 
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or 


Q= 


— t b , () 


T D • 

Io& a 


(D25) 


where <r is the distance between the plates. Equating these 
heat flows results in the following system, of equations: 


- (Tb..- Tb. t)=h. T B ,,)=hi (T Bii -Tii ff>26) 

< 7 


From simultaneous solutions of equations (D21), (D23), and 
(D25) it is found that 



B'rom a simultaneous solution of equation (D24) and equation 
(D24) with the inner boundary condition applied and use of 
equation (D25), it is found that 

i D 

log ‘7T 

T B =T B , i +—£(T B , 0 -T B . i ) (24) 

From equations (D23) and (D25) and from equations (D21) 
and (D25) there are obtained 


A simultaneous solution of equations (D26) gives 


(26) 


(27) 


From equations (26) and (27) it follows that 

Tb..-T bj = (D27) 

1+ ahX {h,+h<) 
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and 


T B . t = 




T b , o 






T B .t=T t + 2k ' B — n (T B , 0 -T B , t ) (22) 

D t h ( log.g 


T B ,,= T vt 2kB - n (T B ',-T B , t ) (23) 

D.h. log.g 


Equations (22) to (25) express the blade temperatures in 
terms of known quantities. 

Section of blade approximated by parallel plates. — The 
derivation of the equations expressing the blade temperatures 
at the inner and outer edges of a blade section approximated 
by a region between parallel plates is given as follows 
(fig. 4(d)): 

(2= heat entering from the hot gas = h 0 { T ti , — T B , 0 ) 

<2= heat leaving to the coo\ant=h ( {T Bit — T t ) 

* k 

(2= heat flowing through the section=— (T Bt0 — T Bit ) 
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